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. Abstract. We present the results of a survey for damped (DLA, log A^(H i) > 20.3) and sub-damped Lyman-a systems 
(19.5 < log A'(H i) < 20.3) at z > 2.55 along the lines-of-sight to 77 quasars with emission redshifts in the range 4 < Zcm < 6.3. 

, Intermediate resolution (R ~ 4300) spectra have been obtained with the Echellette Spectrograph and Imager (ESI) mounted 

^ ■ on the Keck telescope. A total of 100 systems with log A'(H i) > 19.5 are detected of which 40 systems are damped Lyman-o- 

' systems for an absorption length of AX = 378. About half of the lines of sight of this homogeneous survey have never been 

fT^ , investigated for DLAs. We study the evolution with redshift of the cosmological density of the neutral gas and find, consistently 

OO ' with previous studies at similar resolution, that Qdla.Hi decreases at z > 3.5. The overall cosmological evolution of Q.h, shows a 

' peak around this redshift. The H i column density distribution for log A'(H i) > 20.3 is fitted, consistently with previous surveys, 

. with a single power-law of index a — 1.8±0.25. This power-law overpredicts data at the high-end and a second, much steeper, 

' power-law (or a gamma function) is needed. There is a flattening of the function at lower H i column densities with an index of 

^\ ^ a 1.4 for the column density range log ^(H i) = 19.5-21. The fraction of H i mass in sub-DLAs is of the order of 30%. The 

. H I column density distribution does not evolve strongly from z ~ 2.5 to z ~ 4.5. 
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1. Introduction order of log A^(H i) ~ 18, the neutral phase corresponds to 

log A^(H i) > 19.5 (Viegas 1995). The column density defining 

The amount of neutral gas in the Universe is an important n j j j t /t^t an » u u » i 

° ^ the so-called damped Lyman-a (DLA) systems has been taken 

ingredient of galaxy formation scenarios because the neutral . , , ,r.u n ^ on a u tu- j » »u 

° ° to be log A'(H i) > 20.3 because this coiTesponds to the criti- 

phase of the intergalactic medium is the reservoir for star- , f j f ^ f -.■ /m, » i 

^ ° cal mass surface density limit for star formation (Wolfe et al. 

formation activity in the densest places of the universe where ,no<;Nu » i u »u ■ ^ ^ a^u e^u j 

^ ^ 1986) but also because the equivalent width of the correspond- 

galaxies are to be formed. It is therefore very important to make • u • ■ ^ t u e ^^u » • 

° J f ijig absorption is appropriate for a search of these systems in 

a census of the mass in this phase and to determine its cosmo- , , .• ^ n-u t i j u u 

^ low resolution spectra. Therefore several definitions have been 

logical evolution (see e.g. Peroux et al. 2001). j j /-iDla • »u a tu • n,r a 

" V & / introduced. 12^ is the mass density of baryons in DLA sys- 

The gas with highest H I column density is detected through , ^ ■ ,. ., , ,,,,, s ^„ ^ ^h, 

■> tems, defined arbitrarily as systems with log A^(Hi) > 20.3. Qg 



damped Lyman-o' absorptions in the spectra of remote quasars. 
Although damped wings are seen for column densities of the 



is the mass density of neutral hydrogen in all systems : DLAs, 
Lyman limit systems (LLS) and the Lyman-o- forest. The mass 
W oM^rint '^sts to: R. Guimaraes, e-mail: density of H I in the Lyman-o- forest is negligible because the 
rguimara@eso org slope of the H i column density distribution is larger than -2 
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Research in Astronomy, a scientific partnership among the California systems is very difficult to derive directly because the Lyman-o 
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Aeronautics and Space Administration. 



2 



R. Guimaraes et al.: Damped and sub-damped Lyman-cf absorbers in z > 4 QSOs 



However, Qg ' is not easily related to physical quantities as 
the LLS with log < 19.5 are at least partly ionized when the 
ones with log > 19.5 are not (see e.g. Meiring et al. 2008). On 
the contrary, as emphasized by Prochaska et al. (2005), here- 
after PHW05, the mass density of the neutral phase, il™"', is 
a good indicator of the mass available for star-formation and 
should be prefered instead. Note that Q"''"' is not equal to 
The column density limit at which the gas is mostly neutral can- 
not be defined precisely but should lie between log A^(H i) = 19 
and 19.5. In any case, a conservative position is to consider that 
all systems above 19.5 are neutral. 

Whether or not the mass of the neutral gas in the systems 
with 19.5 < log A^(H i) < 20.3 (the so-called sub-DLAs or 
super-LLS) is negligible has been the source of intense discus- 
sions in recent years. Note that these discussions are related to 
the mass in the neutral phase only. Indeed, it is known for long 
(e.g. Petitjean et al. 1993) that the total mass associated with 
the Lyman limit systems is larger than that of DLAs. Indeed the 
gas in the LLS phase is mostly ionized and located in extended 
halos whereas DLAs are located in dense and compact regions. 
Peroux et al. (2003), hereafter PMSI03, have been the first to 
consider the sub-DLAs as an important reservoir of neutral gas. 
They claim that at z > 3.5, DLAs could contain only 50% of 
the neutral gas, the rest being to be found in sub-DLAs. When 
correcting for this, they find that the comoving mass density 
shows no evidence for a decrease above z - 2. PHW05 ques- 
tionned this estimate. They use the Sloan Digital Sky Survey to 
measure the mass density of predominantly neutral gas Qg''"*. 
They find that DLAs contribute >80% of n™"' at all redshift. 
Uncertainties are very large however and the same authors es- 
timate that the systems with log A^(H i) > 19 (the super-LLS) 
could contribute 20-50 % of q!?'. Therefore, the question of 
what is the contribution of super-LLS to O™"' is not settled yet. 

In addition, the evolution of O™"" at the very high redshift, 
z > 4, is not known yet. PHW05 claim that there is no evolution 
of Qg'''^ for z > 3.5 but they caution the reader that results for 
z > 4 should be confirmed with higher resolution data. The 
reason is that the Lyman-a forest is so dense at these redshifts 
that it is very easy to misidentify a strong blend of lines with a 
DLA. Therefore Qg '"^ can be easily overestimated. 

In this paper we present the result of a survey for DLAs 
and sub-DLAs at high redshift (z > 2.55) using intermedi- 
ate resolution data. We identify a total of 100 systems with 
log A^(H i) < 19.5 of which 40 are DLAs over the redshift range 
2.88 < Zahs ^ 4.74 along 77 lines-of-sight towards quasars with 
emission redshift 4 < Zem < 6.3. The sample and data reduction 
are presented in Section 2. In Section 3 we describe the proce- 
dures used to select the absorption systems. Section 4 analyses 
statistical quantitities characterizing the evolution of DLAs and 
sub-DLAs and discusses the cosmological evolution of the neu- 
tral gas mass density. Conclusions are summarized in Section 
5. Throughout the paper, we adopt 0,^ - 0.3, Q.a = 0.7 and Hq 
= 72kms-i. 

2. Observation and data reduction 

Medium resolution (R ~ 4300) spectra of all z > 3 quasars dis- 
covered in the course of the DPOSS survey (Digital Palomar 



Observatory Sky Survey; see, e.g., Kennefick et al. 1995, 
Djorgovski et al. 1999 and the complete listing of QSOs avail- 
able at http://www.astro.caltech.edu/~george/z4.qsos) have 
been obtained with the Echellette Spectrograph and Imager 
(ESI, Sheinis et al 2002) mounted on the KECK II 10 m tele- 
scope. In total, 99 quasars have been observed, 57 of which 
already reported in the literature (see Table [T] for details). We 
provide in Table [T] a summary of the observation log for the 99 
quasars. Columns 1 to 8 give, respectively, the quasar's name, 
the emission redshift, the apparent R magnitude, the J2000 
quasar coordinates, the date of observation, the exposure time 
and the notes. 

The echelle mode allows to cover the full wavelength range 
from 3900 A to 10900 A in ten orders with ~ 300 A overlap 
between two adjacent orders. The instrument has a spectral dis- 
persion of about 11.4 km s"' pixeL' and a pixel size ranging 
from 0.16 A pixel ' in the blue to 0.38 A pixel ' in the red. 
The 1 arcsec wide slit is projected onto 6.5 pixel resulting in a 
R ~4300 spectral resolution. 

Data reduction followed standard procedures using IRAF 
for 70% of the sample and the programme makee for the re- 
maining. For the IRAF reduction, the procedure was as follows. 
The images were overscan corrected for the dual-amplifier 
mode. Each amplifier has a different baseline value and differ- 
ent gain which were corrected by using a script adapted from 
LRIS called esibias. Then all the images were bias subtracted 
and corrected for bad pixels. The images were divided by a nor- 
malized two-dimensional flat-field image to remove individual 
pixel sensitivity variations. The flat-field image was normal- 
ized by fitting its intensity along the dispersion direction using 
a high order polynomial fit, while setting all points outside the 
order aperture to unity. The echelle orders were traced using the 
spectrum of a bright star Cosmic rays were removed from all 
two-dimensional images. For each exposure the quasar spec- 
trum was optimally extracted and background subtracted. The 
task apall in the IRAF package echelle was used to do this. 
The CuAr lamps were individually extracted using the quasar's 
apertures. Lines were identified in the arc lamp spectra by us- 
ing the task ecidentify and a polynomial was fitted to the line 
positions resulting in a dispersion solution with a mean RMS 
of 0.09 A. The dispersion solution computed on the lamp was 
then assigned to the object spectra by using the task dispcor. 
Wavelengths and redshifts were computed in the heliocentric 
restframe. The different orders of the spectra were combined 
using the task scombine in the IRAF package. It is important 
to note that the signal-to-noise ratio drops sharply at the edges 
of the orders. We have therefore carefully controlled this pro- 
cedure to avoid any spurious feature. 

The SNR per pixel was obtained in the regions of the 
Lyman-a forest that are free of absorption and the mean SNR 
value, averaged between the Lyman-a and Lyman-yS QSO emis- 
sion lines, was computed. We used only spectra with mean 
SNR > 10. For simplicity, we excluded from our analysis broad 
absorption line (BAL) quasars. We therefore used 77 lines-of- 
sight out of the 99 available to us. 

The continuum was automatically fitted (Aracil et al. 2004; 
Guimaraes et al. 2007) and the spectra were normalized. We 
checked the normalization for all lines-of-sight and manu- 
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ally corrected for local defects especially in the vicinity of 
the Lyman-Qf emission lines and when the Lyman-o' forest is 
strongly blended. 

Metallicities measured for twelve z > 3 DLAs observed 
along five lines-of-sight of this sample have been published by 
Prochaska et al. (2003b); see also Prochaska et al. (2003a). 



3. Identification of Damped and sub-Damped 
Lyman-Qf systems 

We used an automatic version of the Voigt profile fitting routine 
VPFIT (Carswell et al. 1987) to decompose the Lyman-a forest 
of the spectra in individual components. As usual, we restricted 
our search outside of 3000 km s ' from the QSO emission red- 
shift. This is to avoid contamination of the study by proximate 
effects such as the presence of overdensities around quasars 
(e.g. RolUnde et al. 2005, Guimaraes et al. 2007). 

From this fit we could identify the candidates with 
log A^(H i) - error > 19.5. We then carefully inspected each 
of these candidates individually with special attention to the 
following characteristics of DLA absorption lines: 

- the wavelength range over which the line is going to zero; 

- the presence of damped wings; 

- the identification of associated metal lines when possible. 

The redshift of the H i absorptions was adjusted carefully using 
the associated metal absorption features and the final H i col- 
umn density was then refitted using the high order lines in the 
Lyman series when Lyman series are covered by our spectrum 
and not blended with other lines. After applying the above cri- 
teria, 100 DLAs and sub-DLAs with log N(li i) > 19.5 were 
confirmed. A total of 65 systems were not previously used for 
f^DLA.Hi /^^sub-DLA.Hi estimations, 21 of which are DLAs and 
44 are sub-DLAs. Their characteristics are given in Table |2] 
Columns 1 to 7 give, respectively, the QSO's name; the emis- 
sion redshift estimated as the average of the determinations 
from the peak of the Lya emission and from the peak of a 
Gaussian fitted to the CIV emission line; the minimum red- 
shift along a DLA/sub-DLA could be detected; the maximum 
redshift along a DLA/sub-DLA could be detected; the redshift 
of the DLA/sub-DLA; the DLA/sub-DLA H i column den- 
sity; associated detected metal lines and notes. Lyman series 
and selected associated metal lines are shown in the Appendix 
"List of Figures". The QSO lines of sight of our sample along 
which we detect no DLA and/or sub-DLA are listed in Table [3] 
Comments on DLA/sub-DLA systems differences between our 
measurements and measurements by others are in the Appendix 
"Notes on Individual Systems". 

The metal lines have been searched for using a search 
list of the strongest atomic transitions given in Table H] The 
corresponding absorptions have been fitted using the package 
VPFIT. A full account of this metallic column densities and the 
corresponding abundances is out of the scope of this paper and 
will be presented elsewhere (Guimaraes et al., in preparation). 



4. Analysis 

Using the procedure described in the previous Section, we de- 
tect 100 systems with log (H i) > 19.5, out of which 40 are 
DLAs. We use this sample to investigate the characteristics of 
the neutral phase over the redshift range 2.5 < z < 5. We com- 
pare in Fig[T]the redshift sensitivity function of our survey with 
the redshift sensitivity function of previous surveys computed 
by PMSI03. Although the redshift path of our survey is much 
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Fig. 1. Solid curve shows the redshift sensitivity function of our 
survey for all lines-of-sight, the filled grey histogram is only 
for unpublished lines-of-sight. Dashed curve shows the redshift 
sensitivity function of previous surveys computed by PMSI03. 



smaller than that of the SDSS survey (PHW05), it is at z > 3.5 
similar to the surveys by Peroux et al. (2001, 2003). Note that 
our survey is homogeneous and at spectral resolution twice or 
more larger than previous surveys. The H i column density dis- 
tribution of the 100 DLAs and sub-DLAs measured in this work 
and that of the systems with log A^(H i) > 20.3 in PMSI03 are 
shown in Fig.|2] We are confident that we do not miss a large 
number of sub-DLAs down to the above limit. 

To give a global overview of the survey, we plot in Fig. [3] 
log A^(H i) versus Zabs for the 65 unpublished damped/sub-DLA 
absorption systems. In the same figure we show for comparison 
the data points from the Peroux et al. (2001) survey. 

In Figure H] we plot for comparison, as a function of red- 
shift, the difference between the H i column densities measured 
for the same systems by us and by either Peroux et al. (2005) 
from high-resolution data or Prochaska & Wolfe (2009) from 
SDSS data. The measurements are consistent within errors. 
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Fig. 2. Histogram of the H i column densities measured for the 
100 DLAs and sub-DLAs with log A^(H i) > 19.5 detected in 
our survey (solid-line histogram). The dashed-line histogram 
represents the H i column densities measured by PMSI03 with 
log A^(H i) > 20.3. 



Fig. 4. The H i column density difference between our mea- 
surements and those by either Peroux et al. (2005) (triangles) 
or Prochaska & Wolfe (2009) (squares) in systems common to 
different surveys is plotted versus redshift. 
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Fig. 3. The logarithm of the H i column density measured in our 
unpublished systems (circles) is plotted versus the redshift. The 
data points of Peroux et al. (2001) are shown for comparison 
with crosses. 



4.1. Column Density Distribution Function 

The H I absorption system frequency distribution function is 
defined as: 



f(N,X)dNdX 



AN X AXi 



dNdX 



(1) 



where msys is the number of absorption systems with a col- 
umn density comprized between - AA^/2 and + AN/2 and 
observed over an absorption distance interval of AX. The to- 
tal absorption distance coverage, 2" AXi, is computed over the 
whole sample of n QSO lines-of-sight. The absorption distance, 
X, is defined as 



Xiz) 



(1 +zfE{z)dz 



(2) 



where E(z) - [QmCI + z)^ + Qa]"''^- For one line-of-sight, 
AX(z) = X(zm-dx)-X(Zmm), with z^ax being the emission redshift 
minus 3000 km s ' and Zmin is the redshift of an H i Lyman-a 
line located at the position of the QSO Lyman-/? emission line. 
In Figure |5] we show the function f{N, X) obtained from our 
statistical sample over the redshift range z - 2.55 - 5.03 and 
for log A^(H i) > 19.5. It can be seen that there is no break at 
the low column density end, between log A^(H i) - 19.5 and 
20.6. This makes us confident that we are complete down to 
log A^(H i) = 19.5. The vertical bars indicate Icr errors. The 
horizontal bars indicate the bin sizes plotted at the mean col- 
umn density for each bin. PHW05 results in the redshift range 
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Fig. 5. Frequency distribution function over the redshift range 
z = 2.55 -5.03. The dashed green, soHde black and dotted blue 
lines are, respectively, a power-law, a double power-law and a 
gamma function fits to the data. The solid red and dashed black 
lines are, respectively, a power-law and gamma function fits to 
the data obtained by PHW05. 



z - 2.2 - 5.5 and for logA^(H i) > 20.3 are overplotted in 
the same figure. Although our data points are consistent within 
about Icr with those of PHW05, it seems that the overwhole 
shape of the function is flatter in our data. Note that we do not 
detect any system with log A^(H i) > 21.25. Data points from 
Peroux et al. (2005), hereafter PDDKM05, are also overplot- 
ted. Our point at logA^(H i) - 19.5 is consistent with theirs. 

A power-law, a gamma function and/or a double power- 
law are usually used to fit the frequency distribution. It is ap- 
parent from Fig. |5]that a power-law (of the form f{N,X) = 
K X ") fits the function nicely over the column density range 
19.5 < log A^(H i) < 21. The index of this power spectrum (see 
Table |5]i is larger (a — 1.4) than what is found by PHW05 
but over a smaller column density range log A^(H i) > 20.3. 
The discrepancy is apparently due to the difference in the col- 
umn density ranges considered by both studies. If we restrict 
our fit to the same range as PHW05 we find an index of a — 
1.8 ±0.25 which is consistent with the results of PHW05. We 
note that PDDKM05 akeady mentioned that the small end of 
the column density distribution is flatter than a - -2. 

There is a large deficit of high column density systems in 
our survey compared to what would be expected from the sin- 
gle power-law fit. This has been noted before and discussed 
in detail by PHW05. A double power-law was used to fit our 
full sample with better result (see Table|5]l. However, the sharp 
break in the function at log A^(H i) ~ 21 suggests a gamma func- 



tion of file form, f{N, X) = Kx xe^ (Pei & Fafl 1995) 

should better describe the data. 

We have calculated the frequency distribution function, 
f(N,X), in different redshift bins of equal distance path : 
2.55-3.40; 3.40-3.83 and 3.83-5.03. Results are shown in Fig. 
|6l The functions are fitted as described above and fit results are 
given in Table |5] We find that the function do not evolve much 
with redshift. This is consistent with the finding by PHW05 
that the global characteristics of the function are not changing 
much with time. There is however a tendancy for a flattening of 
the function which may indicate that the number of sub-DLAs 
relatively to other systems is larger at lower redshift. Although 
we do not think this is the case because we have used a con- 
servative approach, part of this evolution at the highest redshift 
could possibly be a consequence of loosing the sub-DLAs in 
the strongly blended Lyman-a forest at z > 4. We note also a 
slight decrease of the number of systems with log A^(H i) > 20.5 
at the highest redshifts. This is consistent with the finding by 
PDDKM05 that the relative number of high column density 
DLAs decreases with redshift. 

Another way to look at these variations is to compute the 
redshift evolution of the number of (sub)DLAs per unit path 
length. The observed density of systems is defined as 



(sub)DLA,s 



f(N,X)dN 



(3) 



Results obtained during the present survey together with 
those of PWH05 and PDDKM05 are plotted in Figure |7] The 
important feature of this plot is that the number density of 
DLAs peaks at z ~ 3.5. In addition, the ratio of the number 
of sub-DLAs to the number of DLAs is larger at redshift <3.5 
compared to higher redshifts. 

4.2. Neutral hydrogen cosmological mass density, Qhi 

The comoving mass density of neutral gas is given by 
Hofimn 2;,M(Hi) 



c po AX(z) 



(4) 



where the density is in units of the current critical density 
Po; niu is the mass of the hydrogen atom, fj. = 1.3 is the mean 
molecular weight, AX the absorption length and the summa- 
tion of column densities is done over all absorption systems 
detected in the survey. Results are shown in Fig. [8] and summa- 
rized in Table |6] The three redshift bins considered are defined 
so that the absorption length is equal in each bin. The different 
columns of Table |6] give, respectively, the redshift range, the 
mean redshift, the number of lines-of-sight involved, the num- 
ber of DLAs and sub-DLAs detected over this redshift range, 
the absorption length (calculated using Zmin and Zmax as defined 
in Table|2]l, the resulting H i cosmological densities for, respec- 
tively, DLAs only or both DLAs and sub-DLAs. 

Results from PHW05 and PDDKM05 ai-e also plotted on 
Figure[8] It can be seen that we confirm the decrease of Qhi for 
z > 3 that was noticed by PDDKM05. The measurement from 
the SDSS in this redshift range is higher However, our survey 
is of higher spectral resolution and should in principle be more 
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Fig. 6. H I frequency distribution, /h,(A^,X), in three redshifts bins: 2.55-3.40 (left), 3.40-3.83 (center), and 3.83-5.03 (right). 
Straight Hnes show best fits of a power-law function to the binned data. The dashed curve is the same for a gamma function 
fit. and the dot-dashed curve for a double power-law fit. 



reUable in this redshift range. It seems that the evolution of Qhi 
is a steep increase from z = 2 to z = 3 and then a slightly flatter 
decrease up to z = 5. The inclusion of the sub-DLAs does not 
change this picture as sub-DLAs contribute to a maximum of 
about 30% to the total H i mass. The contribution by sub-DLAs 
is better seen in Fig. |9] where we plot the cumulative density 
versus the maximum H i column density considered. As noted 
already by numerous authors, the discrepancy of measurements 
at z < L5 is still a problem. 

It can be also seen from the Figure [8] that Qh, are lower 
than Qsteiiar, the mass density in stars in local galaxies. We find 
for the ratio of the peak value of fie, to i^steiiai for this work R 
~ 0.45. Previous surveys, PDDKM05 and PHW05, have found 
for the ratio 0.37 and 0.40 respectively. 

5. Discussion 

We have presented the results of a survey for damped and sub- 
damped Lyman-ff systems (logA^(H i) > 19.5) at z > 2.55 
along the lines-of-sight to 77 quasars with emission redshifts 
in the range 4 > Zem ^ 6.3. In total 99 quasars were observed 
but 22 lines-of-sight were not used because of not enough 
SNR and/or because of the presence of broad absorption lines. 



Intermediate resolution (R ~ 4300) spectra have been obtained 
with the Echellette Spectrograph and Imager (ESI) mounted 
on the Keck telescope. The damped Lyman-a absorptions are 
identified on the basis of (i) the width of the saturated absorp- 
tion, (ii) the presence of damped wings and (iii) the presence 
of metals at the corresponding redshift. The detection is run 
automatically but all lines are verified visually. A total of 100 
systems with log A^(H i) > 19.5 are detected of which 40 sys- 
tems are Damped Lyman-a systems (log A^(H i) > 20.3) for an 
absorption length of AX - 378. Spectra are shown in Appendix. 

PHW05 have derived from SDSS data that the cosmolog- 
ical density of the neutral gas increases strongly by a factor 
close to two from z ~ 2 to z ~ 3.5. Beyond this redshift, 
measurements are more difficult because the Lyman-a forest is 
dense. Our measurements should be more reliable because of 
better spectral resolution. We show, consistently with the find- 
ings of PDDKM05, that the cosmological density of the neutral 
gas decreases at z > 3.5. The overall cosmological evolution 
seems therefore to have a peak at this redshift. 

We find that the H i column density distribution does not 
evolve strongly from z ~ 2.5 to z ~ 4. The one power-law fit in 
the range log A^(H i) > 20.3 gives an index of a = -L80+0.25, 
consistent with previous determinations. However, we find that 
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Fig. 7. Number density of absorbers vs. redshift. Red circles 
and inverse triangles are from this work for DLAs and sub- 
DLAs, respectively. The values obtained by PWH05 (blue in- 
verse triangles) and PDDKM05 (green squares) for DLAs are 
overplotted. The green diamonds are the values obtained for 
log N(H i) > 19.0 by PDDKM05. 



the fit over the column density range log A^(H i) = 19.5-21 is 
quite flat (a ~ 1.4). This probably indicates that the slope at 
the low-end is much flatter than -2. This power-law overpre- 
dicts data at the high-end and a second much steeper power-law 
(or a gamma function) is needed. The fraction of H i mass in 
sub-DLAs is of the order of 30%. Our data do not support the 
claim by PDDKM05 that the incidence of low column-density 
sytems is larger at high redshift. The number density of sub- 
DLAs seems to peak at z ~ 3.5 as well. 

It is apparent that statistical errors are still large in our sur- 
vey. It would be therefore of first importance to enlarge the 
sample of (sub)DLAs at high redshift. For this we need to go to 
fainter quasars however. The advent of X-shooter, a new gen- 
eration spectrograph at the VLT with a spectral resolution of 
R - 6700 in the optical, should allow this to be done in a rea- 
sonable amount of observing time. 
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Fig. 8. Cosmological evolution of the H i mass density. For 
DLAs: red squares are the results of this work, green squares 
are from PDDKM05, blue inverse triangles from PHW05 and 
black triangles from Rao et al. (2005). For sub-DLAs: red cir- 
cles are this work (for systems with column densities > 19.5), 
green circles obtained from PDDKM05 (for systems with col- 
umn densities > 19.0) and blue circles obtained from PHW05 
(as determined from the single power-law fit to the LLS fre- 
quency distribution function). The black circle are the mass 
density in stars in local galaxies (Fukugita et al. 1998). 
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QSO 




R 


RA(2000) 


DEC(2000) 


Obs. Date 


Exp. Time (min) 


Note* 


PSS0007-I-2417 


4.050 


18.69 


00 07 38.7 


+24 17 25.0 


2000 Sep 06 


60 




PSS0014+3032 


4.470 


18.81 


00 14 42.8 


+30 32 03.0 


2000 Sep 06 


60 




PSS0052+2405 


4.280 


18.30 


00 52 06.8 


+24 05 39.0 


1999 Dec 30 


40 




PSSOl 17+1552 


4.244 


18.60 


01 17 31.2 


+15 52 16.4 


2000 Sep 04 


50 


2 


PSS0118+0320 


4.232 


18.50 


01 18 52.8 


+03 20 50.0 


2000 Sep 06 


60 


3 


PSS0121+0347 


4.127 


17.86 


01 21 26.2 


+03 47 07.0 


1999 Dec 30 


45 


3 


SDSS0127-0045 


4.084 


18.37 


01 27 00.7 


-00 45 59.4 


2001 Jan 02 


60 




PSS0131+0633 


4.430 


18.24 


01 31 12.2 


+06 33 40.0 


2000 Sep 04 


60 


1,2 


PSS0133+0400 


4.154 


17.86 


01 33 40.4 


+04 00 59.0 


1999 Dec 29 


40 


1,2,3 


PSS0134+3307 


4.536 


18.82 


01 34 21.6 


+33 07 56.5 


2000 Sep 06 


60 


1,2 


PSS0207+0940 


4.136 


18.63 


02 07 03.5 


+09 40 59.0 


2000 Sep 04 


75 




SDSS0206+1216 


4.810 


21.51 


02 06 51.4 


+ 12 16 24.4 


2002 Dec 07 


30 




PSS0209+0517 


4.194 


17.36 


02 09 44.7 


+05 17 14.0 


1999 Dec 29 


35 


1,2,3 


SDSS0210-0018 


4.700 


20.74 


02 10 43.2 


-00 18 18.5 


2002 Dec 07 


40 


4,5 


PSS0211+1107 


3.975 


18.12 


02 11 20.1 


+ 11 07 16.0 


1999 Dec 29 


70 




SDSS021 1-0009 


4.900 


22.04 


02 11 02.7 


-00 09 10.3 


2002 Dec 07 


20 


1,2 


SDSS023 1-0728 


5.410 


21.54 


02 31 37.6 


-07 28 54.5 


2002 Dec 06 


90 


4,5 


PSS0244-0108 


3.990 


19.00 


02 44 57.2 


-01 08 08.7 


2000 Sep 06 


40 


4,5 


PSS0248+1802 


4.422 


18.40 


02 48 54.3 


+18 02 50.3 


1999 Dec 30 


60 


1,2 


PSS0320+0208 


3.960 


18.74 


03 20 42.7 


+02 08 16.0 


1999 Dec 30 


55 




SDSS0338+0021 


5.020 


21.68 


03 38 29.3 


+00 21 56.5 


2002 Dec 06 


90 


1, 2, 4, 5 


SDSS0338-RD657 


4.960 


23.00 


03 38 31.3 


+00 18 07.7 


2001 Jan 02 


120 




PSS0452+0355 


4.420 


18.80 


04 52 51.5 


+03 55 58.0 


1999 Dec 30 


60 




PSS0747+4434 


4.435 


18.06 


07 47 50.0 


+44 34 16.0 


1999 Dec 30 


90 


1,2,5 


SDSS0756+4104 


5.090 


21.70 


07 56 18.0 


+41 04 10.6 


2001 Mar 26 


150 


4 


PSS0808+5215 


4.510 


18.82 


08 08 49.5 


+52 15 16.0 


2000 Apr 28 


70 


5 


SDSS08 10+4603 


4.074 


18.67 


08 10 54.7 


+46 03 55.2 


2001 Mar 24 


60 




PSS0852+5045 


4.200 


19.00 


08 52 27.4 


+50 45 11.0 


2000 May 13 


60 




PSS0926+3055 


4.190 


17.31 


09 26 36.3 


+30 55 06.0 


1999 Dec 29 


40 




SDSS094 1+5947 


4.820 


20.66 


09 41 08.4 


+59 47 25.8 


2002 Dec 06 


40 


4,5 


BR0945-0411 


4.130 


18.80 


09 47 49.6 


-04 25 15.1 


2000 May 14 


60 




PSS0950+5801 


3.973 


17.38 


09 50 14.0 


+58 01 38.0 


1999 Dec 29 


40 


4,5 


BR095 1-0450 


4.350 


18.90 


09 53 55.7 


-05 04 19.5 


2000 May 15 


53 


2 


PSS0955+5940 


4.340 


17.84 


09 55 11.3 


+59 40 32.0 


1999 Dec 29 


56 


4,5 


PSS0957+3308 


4.283 


17.59 


09 57 44.5 


+33 08 23.0 


1999 Dec 29 


55 


5 


BRI1013+0035 


4.380 


18.80 


10 15 49.0 


+00 20 19.0 


2000 May 15 


60 


2 


PSS1026+3828 


4.180 


18.93 


10 26 56.7 


+38 28 43.0 


2000 Apr 29 


95 


5 


BR1033-0327 


4.509 


18.50 


10 36 23.7 


-03 43 20.0 


2000 May 13 


20 




PSS 1039+3445 


4.390 


19.20 


10 39 19.3 


+34 45 10.9 


2001 Apr 18 


150 


5 


SDSS1044-0125 


5.740 


25.10 


10 44 33.0 


-01 25 03.1 


2000 Apr 28 & 2001 Jan 01,02 


330 




SDSS 1048+4637 


6.230 


22.40 


10 48 45.0 


+46 37 18.3 


2003 Jun 02 


90 


5 


PSS 1048+4407 


4.450 


19.50 


10 48 46.6 


+44 07 12.7 


2000 Apr 29 


20 




PSS1057+4555 


4.126 


17.70 


10 57 56.4 


+45 55 51.9 


1999 Dec 30 


40 


1,2,5 


PSS1058+1245 


4.330 


18.00 


10 58 58.5 


+ 12 45 55.0 


1999 Dec 30 


75 




PSSl 118+3702 


4.030 


18.76 


11 18 56.2 


+37 02 53.9 


2001 Mar 24 


60 


5 


PSS 1140+6205 


4.509 


18.73 


1140 09.6 


+62 05 23.3 


2000 May 14 


60 


4,5 


PSS1159+1337 


4.081 


18.50 


11 59 06.5 


+ 13 37 37.8 


2000 Feb 10 


55 


1,2,5 


SDSS1204-0021 


5.030 


20.82 


12 04 41.7 


-00 2149.6 


2003 Jun 04 


60 


4,5 


SDSS1208+0010 


5.273 


22.75 


12 08 23.9 


+00 10 28.9 


2001 Mar 24 


120 




PSS 1226+0950 


4.340 


18.78 


12 26 23.8 


+09 50 03.7 


2001 Mar 24 


80 




PSS 1247+3406 


4.897 


20.40 


12 49 42.2 


+33 49 54.0 


2001 Mar 26 


60 




PSS1248+3110 


4.346 


18.90 


12 48 20.2 


+31 10 44.0 


2000 May 13 


60 




PSS 1253-0228 


4.007 


19.40 


12 53 36.3 


-02 28 08.0 


2000 Apr 29 


55 


1,2 


SDSS 13 10-0055 


4.151 


18.85 


13 10 52.6 


-00 55 31.8 


2001 Mar 24 


50 


4,5 
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Table 1. continued. 



QSO 




R 


RA(2000) 


DEC(2000) 


Obs. Date 


Exp. Time (min) 


Note* 


PSS 13 15+2924 


4.180 


18.48 


13 15 39.6 


+29 24 39.8 


2001 Mar 24 


85 




PSS1317+3531 


4.369 


19.10 


13 17 43.3 


+35 31 33.1 


2001 Mar 26 


55 


2 


J1325+1123 


4.400 


18.77 


13 25 12.6 


+ 11 23 32.8 


2001 Apr 18 


120 


5 


PSS 1326+0743 


4.123 


17.30 


13 26 11.9 


+07 43 59.0 


2000 Dec 29 


60 




PSS1339+5154 


4.080 


18.70 


13 39 13.0 


+51 54 04.0 


2000 May 15 


60 


5 


PSS 1347+4956 


4.560 


17.90 


13 47 43.4 


+49 56 21.0 


2000 Feb 10 


50 


5 


PSS1401+4111 


4.026 


18.62 


14 01 32.8 


+41 11 49.9 


2000 Apr 28 


46 


4, 5 


PSS 1403+4126 


3.862 


18.92 


14 03 55.7 


+41 26 16.2 


2000 May 15 


90 


4, 5 


PSS1418+4449 


4.280 


18.40 


14 18 31.8 


+44 49 37.0 


2000 Feb 29 


60 


4,5 


PSS 1430+2828 


4.306 


19.30 


14 30 31.9 


+28 28 34.1 


2001 Mar 26 


55 


2 


PSS 1432+3940 


4.292 


18.60 


14 32 24.9 


+39 40 24.0 


2000 May 13 


60 


4, 5 


PSS 1435+3057 


4.350 


19.30 


14 35 23.3 


+30 57 16.3 


2001 Apr 18 


70 


2 


PSS 1443 +2724 


4.406 


19.30 


14 43 31.2 


+27 24 37.0 


2000 Apr 29 


30 


2 


PSS1443+5856 


4.270 


17.80 


14 43 40.8 


+58 56 53.0 


2000 Apr 28 


40 


4,5 


PSS1458+6813 


4.291 


18.67 


14 58 31.7 


+68 13 05.2 


2000 Apr 28 


60 




PSS1500+5829 


4.224 


18.60 


15 00 07.7 


+58 29 38.0 


2000 May 14 


60 




PSS 1506+5220 


4.180 


18.10 


15 06 54.6 


+52 20 05.0 


2000 Feb 10 


50 


4,5 


GB1508+5714 


4.304 


18.90 


15 10 02.2 


+57 03 04.9 


2000 May 13 


75 


2 


PSS1531+4157 


4.200 


19.00 


15 31 29.4 


+45 17 07.9 


2001 Apr 18 


128 




PSS 1535+2943 


3.972 


18.90 


15 35 53.9 


+29 43 13.0 


2000 May 1 5 & 2000 Sep 05 120 




PSS1543+3417 


4.390 


18.40 


15 43 40.4 


+34 17 45.0 


2000 Apr 30 


76 


5 


PSS1554+1835 


3.990 


18.90 


15 35 53.9 


+29 43 13.0 


2000 May 15 & 2000 Sep 04 120 




PSS 1555+2003 


4.228 


18.90 


15 55 02.6 


+20 03 25.0 


2000 May 14 


90 




PSS1615+1803 


4.010 


18.42 


16 15 22.9 


+ 18 03 56.4 


2000 Apr 29 


90 




SDSS1630+4012 


6.050 


20.42 


16 30 33.9 


+40 12 09.6 


2003 Jun 02,03,04 


295 




PSS1633+1411 


4.349 


19.00 


16 33 19.7 


+ 14 11 42.6 


2000 Apr 30 


60 


1, 2 


PSS 1646+55 14 


4.084 


18.11 


16 46 56.3 


+55 14 46.7 


2000 Apr 28 


60 


1,2 


VLA1713+4218 


4.230 


19.00 


17 13 56.2 


+42 18 08.6 


2001 Apr 17 


90 




PSS1715+3809 


4.520 


18.56 


17 15 39.5 


+38 09 06.6 


2000 Sep 03 


60 




PSS1721+3256 


4.040 


19.23 


17 21 06.7 


+32 56 35.8 


2000 May 01 


60 


1, 2 


PSS 1723 +2243 


4.514 


18.17 


17 23 23.1 


+22 43 56.4 


2001 Mar 24 


40 




SDSS1737+5828 


4.940 


20.93 


17 37 44.0 


+58 28 25.4 


2001 Mar 24 


90 


4,5 


PSS 1745+6846 


4.130 


19.12 


17 45 50.1 


+68 46 21.0 


2000 May 01 


120 




PSS1802+5616 


4.158 


19.19 


18 02 48.9 


+56 16 51.0 


2000 Sep 05 


90 


1, 2 


PSS2122-0014 


4.114 


19.13 


21 22 07.5 


-00 14 45.0 


2000 Sep 06 


60 


1,2, 4,5 


PSS2154+0335 


4.360 


1 O /I 1 

18.41 


21 54 06.9 


+03 35 40.0 


2000 May 14 


60 


1, 2, 3 


PSS2155+1358 


4.256 


18.50 


21 55 02.1 


+ 13 58 26.0 


2000 May 15 


40 


1,2,3 


PSS2203+1824 


4.375 


18.74 


22 03 43.4 


+ 18 28 14.0 


2000 Sep 04 


60 




PSS2238+2603 


4.031 


18.85 


22 38 41.1 


+26 03 46.0 


2000 Sep 06 


60 




PSS2241+1352 


4.441 


18.69 


22 41 47.7 


+ 13 52 05.9 


2000 Sep 03 


60 


1, 2, 4, 5 


PSS2244+1005 


4.040 


18.92 


22 44 05.5 


+ 10 47 38.0 


2000 Sep 04 


80 




PSS2315+0921 


4.412 


18.96 


23 15 59.2 


+09 21 44.0 


2000 Sep 05 


90 




PSS2322+1944 


4.170 


18.29 


23 22 07.2 


+ 19 44 23.0 


2000 May 1 3 & 2000 Sep 04 111 




PSS2323+2758 


4.180 


18.51 


23 23 41.0 


+27 58 01.0 


2000 Sep 04 


86 




PSS2344+0342 


4.340 


17.87 


23 44 03.2 


+03 42 26.0 


1999 Dec 30 


45 


1,2 



* Quasars from this survey were previously used for Qdla.hi estimations by: (1) Peroux et al. (2001); 

(2) Peroux et al. (2003); (3) Peroux et al. (2005); (4) Prochaska et al. (2005); (5) Prochaska & Wolfe (2009) 
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QSO 




^min 


Zmax 


DLA-subDLA 

Zabs log Nin l) 


Metal Lines 


Note • 


PSS0007+2417 


4.05 


2.758 


4.000 


3.194 


19.50+0.07 


Si II, C IV 












3.497 


20.95+0.10 


Si n, Fe n, Al n 












3.706 


20.65±0.15 


Si n, O I, C n, Si iv, C iv, Fe n, Al n 












3.839 


20.65±0.10 


Si n, O I, C n, C iv, Fe n, Al n 




PSS0117+1552 


4.244 


2.784 


4.192 


3.175 


19.50+0.10 


Civ 




PSSOl 18+0320 


4.232 


2.901 


4.180 


4.128 


19.95±0.15 


Si n, O I, C n, Si rv, Crv, Al n 


3 


PSS0121+0347 


4.127 


2.789 


4.076 


2.977 


19.50+0.15 


no metals " 


3 


SDSSO 127-0045 


4.084 


2.546 


4.033 


3.728 


20.85+0.10 


Si n, O I, C 11, Si iv, C iv, Fe n, Al n 




PSSOl 3 1+0633 


4.430 


3.008 


4.376 


3.173 


19.95±0.10 


Fe n, Al n 


1, 2 










3.689 


19.50+0.10 


C IV 


1, 2 


PSS0133+0400 


4.154 


2.930 


4.102 


3.690 


20.40+0.15 


Siu 


1,2,3 










3.771 


20.45+0.15 


Si n, Si rv, C rv, Fe n, Al n 


1, 2, 3 










3.995 


20.15+0.10 


O I, C II, Si IV, C IV, Fe n 


3 


PSS0134+3307 


4.536 


2.976 


4.481 


3.760 


20.60+0.10 


Si n. Si IV, C iv, Al ii 


1, 2 


PSS0209+0517 


4.194 


2.943 


4.142 


3.666 


20.50+0.10 


Si n, Fe n, Al n 


1, 2, 3 










3.862 


20.30+0.15 


O I, C II, Si II, Al II 


1, 2, 3 


PSS0211+1107 


3.975 


2.847 


3.925 


3.139 


20.05+0.25 


Si n, C IV, Fe ii, Al ii 












3.502 


20.10+0.20 


Si n, C rv, Fe n, Al n 




PSS0747+4434 


4.435 


3.107 


4.381 


3.762 


20.00+0.20 


no metals 


1, 2 










4.019 


21.10+0.15 


O I, C II, Si II, Fe 11, Al ii 


1,2,5 


SDSS0756+4104 


5.09 


3.551 


5.029 


4.360 


20.15+0.10 


no metals 




PSS0808+5215 


4.510 


3.181 


4.455 


3.114 


20.60+0.10 


no metals " 


5 


SDSS0810+4603 


4.074 


2.811 


4.023 


2.956 


20.55+0.25 


Aim 












3.472 


19.90+0.07 


Fe n, Al n 




PSS0950+5801 


3.973 


2.759 


3.923 


3.264 


20.50+0.10 


Si III , Si II, C IV, Fe n, Al n 


4,5 


BR095 1-0450 


4.35 


2.866 


4.296 


3.235 


20.25+0.10 


C IV, Fe II, Al II 












3.856 


20.70+0.20 


Si n. Si IV, C rv, Al n 


2 










4.202 


20.35+0.15 


Si II, O I, C IV 


2 


PSS0955+5940 


4.34 


3.070 


4.287 


3.542 


20.30+0.10 


Si II, C IV, Fe II, Al ii 


4,5 










3.843 


20.00+0.15 


Si IV, Si II, C IV, Al II 


4,5 










4.044 


20.10+0.15 


O I, C II, Si II, Si IV, C IV 




PSS0957+3308 


4.283 


2.969 


4.230 


3.043 


19.65+0.15 


C II, Al n 












3.280 


20.50+0.10 


Fe II, Al n 












3.364 


19.70+0.15 


Al II 












3.900 


19.50+0.10 


no metals 












4.179 


20.60+0.10 


O I, C II, Si II, Si IV, C IV 




BRI1013+0035 


4.38 


3.123 


4.326 


3.738 


20.15+0.20 


no metals 




PSS 1026+3828 


4.18 


2.756 


4.128 


3.339 


19.72+0.08 


no metals over available coverage 












3.865 


19.90+0.10 


C II 




PSS 1039+3445 


4.39 


2.881 


4.336 


4.098 


19.65+0.10 


O I, C II, Si iV 




PSS 1057+4555 


4.126 


2.776 


4.075 


2.909 


20.05+0.10 


Fe II, Al II 


1,2 










3.058 


19.80+0.15 


C iV, Fe II, Al II 


1, 2 










3.164 


19.50+0.20 


Si n, C iV, Fe n, Al n 












3.317 


20.15+0.10 


Si n, Fe n, Al n 


1, 2 


PSS1058+1245 


4.330 


3.110 


4.277 


3.430 


20.40+0.15 


no metals over available coverage 




PSSl 118+3702 


4.03 


2.893 


3.980 


3.698 


19.87+0.08 


O I, C 11, Si 11, Al 11 




PSS1159+1337 


4.081 


2.823 


4.030 


3.726 


20.00+0.10 


C n. Si n. Si iV, C iV 


1,2 


PSS 1247+3406 


4.897 


3.314 


4.838 


4.572 


19.80+0.10 


Si iV, C iV 




PSS 1248+3 110 


4.346 


3.088 


4.293 


3.697 


20.35+0.10 


Si II, Fe II, Al n 












4.075 


20.20+0.15 


O I, C II, Si IV 




PSS 1253-0228 


4.007 


2.868 


3.957 


3.608 


19.95+0.10 


C II, Si IV, C iV, Fe ii 


1,2 


PSS1317+3531 


4.369 


2.998 


4.315 


3.461 


19.90+0.05 


Si n, C iV, Fe n, Al n 




J1325+1123 


4.400 


3.017 


4.346 


3.723 


19.50+0.20 


CiV 
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QSO 


^em 


^min 


^max 


DLA-subDLA 
z,bs log A^(H I) 


Metal Lines 




Note * 










4.133 


19.50±0.20 


no metals 






PSS 1326+0743 


4.123 


2.866 


4.072 


2.919 


19.95+0.10 


Alii 














3.425 


19.90+0.15 


Aln 






PSS 1432+3940 


4.292 


3.105 


4.239 


3.274 


21.00+0.15 


C IV, Fe n, Al n 






PSS 1435+3057 


4.35 


2.942 


4.297 


3.267 


20.05+0.10 


C IV, Al II 














3.516 


20.20+0.10 


Si II, C IV, Al II 














3.778 


19.85+0.10 


C IV, Al II 






PSS 1443+2724 


4.406 


3.095 


4.352 


4.223 


20.90+0.15 


Si IV, Si II, C IV, Fe ii, Al ii 




2 


PSS 1500+5829 


4.224 


3.031 


4.172 


3.585 


19.55+0.10 


no metals over available covera^ 


,g c 












3.915 


20.50+0.10 


Oi 






PSS 1506+5220 


4.18 


3.275 


4.128 


3.223 


20.70+0.10 


Si II, C IV, Fe 11, Al ii, Al iii 






PSS1531+4157 


4.20 


2.868 


4.148 


3.657 


19.55+0.15 


Si II, C IV, Al II 






PSS 1535+2943 


3.972 


2.797 


3.922 


3.202 


20.60+0.10 


Si II, C IV, Fe II, Al ii 














3.762 


20.55+0.15 


O I, C II, Si IV, C IV, Al II 






PSS1554+1835 


3.99 


2.769 


3.940 


2.919 


20.00+0.15 


Si II, C IV, Al II, Al III 






PSS 1555+2003 


4.228 


2.824 


4.176 


3.427 


19.90+0.10 


no metals over available coveraj 


re' 




PSS1633+1411 


4.349 


2.836 


4.296 


2.880 


20.30+0.15 


Fe II, Si II 






PSS 1646+55 14 


4.084 


2.800 


4.033 


2.932 


19.50+0.10 


no metals over available coveraj 


jg c 












4.029 


19.80+0.15 


no metals over available coveraj 


jg c 




PSS1715+3809 


4.52 


3.004 


4.465 


3.341 


21.05+0.20 


Si II, Fe II, Al II 






PSS1723+2243 


4.515 


3.062 


4.460 


3.697 


20.35+0.10 


Si II, C IV, Fe II, Al ii 






SDSS1737+5828 


4.94 


3.383 


4.881 


4.152 


19.85+0.15 


Si II, Al II 














4.741 


20.65+0.15 


Si II, O I, C II 




4,5 


PSS 1745 +6846 


4.13 


2.860 


4.079 


3.706 


19.73+0.09 


Si II, C IV, Al III 






PSS1802+5616 


4.158 


2.821 


4.106 


3.391 


20.25+0.10 


Si IV, Si II, C IV, Fe n, Al ii 




1,2 










3.554 


20.40+0.10 


O I, Si II, C IV, Fe II, Al ii 




1, 2 










3.762 


20.65+0.15 


C II, Si II, Fe II, Al II 




1,2 










3.811 


20.35+0.15 


C II, Si IV, Si II, Al II 




1, 2 


PSS2122-0014 


4.114 


2.903 


4.063 


3.207 


20.20+0.10 


Si II, C IV, Fe II, Al ii 




1, 2, 4, 5 










3.264 


19.90+0.15 


C IV, Al 11 














4.001 


20.15+0.15 


C II, Si IV, C IV 




1,2 


PSS2155+1358 


4.256 


3.064 


4.203 


3.143 


19.90+0.20 


Fe II, Al II 














3.318 


20.75+0.20 


Si II, Fe II, Al II 




1,2 










4.211 


20.00+0.10 


Si II, O I, C II 






PSS2203+1824 


4.375 


2.850 


4.321 


3.610 


19.70+0.15 


Si II, C IV, Al II 














4.107 


19.80+0.15 


C II, Si II, Fe II, Al II 






PSS2238+2603 


4.031 


2.816 


3.981 


3.857 


20.45+0.15 


Si II, O I, C II, Si IV, C IV, Fe ii, Al n 




PSS2241+1352 


A A A t 

4.441 


3.106 


4.387 


3.656 


20.15+0.20 


Al II, Si II 




1 o 
1, 2 










4.281 


20.75+0.15 


Si II, O I, C II, Fe II, Al ii 




1,2, 4,5 


PSS2315+0921 


4.412 


2.863 


4.358 


3.220 


21.25+0.15 


C II, Si II, Fe II, Al II 














3.425 


21.00+0.20 


Si II, C IV, Fe II, Al n 






PSS2322+1944 


4.17 


2.754 


4.118 


2.888 


19.95+0.10 


Fe II, Al II 














2.975 


19.80+0.10 


Fe II, Al II 






PSS2323+2758 


4.18 


2.823 


4.128 


2.952 


19.80+0.10 


no metals 














3.684 


20.60+0.15 


Si II, Fe II 






PSS2344+0342 


4.340 


2.939 


4.287 


3.220 


21.20+0.10 


no metals over available covera^ 


,g c,d 


1,2 










3.884 


19.80+0.10 


no metals over available covera^ 


;e 





" Metals were detected by Peroux et al. (2005) 

* Metals were detected by Prochaska & Wolfe (2009) 

We do not cover the red part of the spectrum A > 6400 A 
Metals were detected by Peroux et al. (2001) 
" Metals were detected by Dessauges-Zavadsky et al. (2003) 

* The same references that appears in Table [T] 
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QSO 




7/. 


z,, 








Nolc 


PSS0014+3032 


4, 


.470 


2. 


.866 


4, 


.415 




SDSS0210-0018 


4, 


.700 


3, 


.177 


4, 


.643 


4,5 


PSS0248+1802 


4, 


.430 


3. 


.118 


4, 


.376 


2 


PSS0452-I-0355 


4, 


.395 


3. 


.115 


4, 


.341 




PSS0852+5045 


4, 


.216 


2, 


.939 


4, 


.164 




PSS0926+3055 


4, 


.198 


2. 


.951 


4, 


.146 




SDSS0941+5947 


4, 


.820 


3, 


.322 


4, 


.762 


4,5 


PSS 1140+6205 


4, 


.509 


3, 


.113 


4, 


.454 


4,5 


SDSS1310-0055 


4, 


.152 


2, 


.830 


4, 


.101 


4, 5 


PSS1339+5154 


4, 


.080 


2. 


.783 


4, 


.029 


5 


PSS1401+4111 


4, 


.026 


2. 


.866 


3, 


.976 


4,5 


PSS 1403+4126 


3, 


.862 


2, 


.866 


3, 


.813 


4,5 


PSS 141 8+4449 


4, 


.323 


2, 


.974 


4, 


.270 


4,5 


PSS 1430+2828 


4, 


.306 


2. 


.811 


4, 


.253 


2 


PSS1458+6813 


4, 


.291 


2. 


.990 


4, 


.238 




GB1508+5714 


4, 


.304 


3. 


.217 


4, 


.251 


2 


PSS 1543+3417 


4, 


.407 


3, 


.071 


4, 


.353 


5 


PSS1615+1803 


4, 


.010 


2, 


.783 


3, 


.960 




PSS1721+3256 


4, 


.040 


2. 


.802 


3, 


.990 


1, 2 


PSS2154+0335 


4, 


.359 


3, 


.400 


4, 


.305 


1,2,3 


PSS2244+1005 


4, 


.040 


2. 


.810 


3, 


.990 





* Quasars from this survey were previously used for nDLA.Hi 

estimations by: (1) Peroux et al. (2001); (2) Peroux et al. (2003) 

(3)Peroux et al. (2005); (4) Prochaska et al. (2005); Prochaska & Wolfe (2009) 



Table 4. Principal absorption metal lines most frequently detected associated to high column density absorption systems. 



Ion 


(A) 


f 


logAof + log[NIN{H)]<i> + 12.00 


Sim 


1206.500 


0.221 


10.87 


Hi 


1215,6701 


0.4162 


14.70 


Nv 


1238.821 


0.152 


10.24 


Nv 


1242.804 


0.0757 


9.93 


Sin 


1260.4223 


0.959 


10.65 


Oi 


1302,1685 


0.0486 


10.67 


Sin 


1304,3702 


0.147 


9.85 


Cn 


1334,5323 


0.118 


10.85 


Sirv 


1393.76018 


0.528 


10.44 


Si IV 


1402.770 


0.262 


10.13 


Sin 


1526,70698 


0.23 


10.11 


Civ 


1548.2041 


0.194 


11.13 


Civ 


1550.7812 


0.097 


10.83 


Fe n 


1608.45085 


0.062 


9.52 


Aln 


1670,7886 


1.88 


9.99 


Aim 


1854.7164 


0.539 


9.49 


Aim 


1862.7895 


0.268 


9.19 


Fe n 


2344.2139 


0.108 


9.92 


Fe n 


2374.4162 


0.0395 


9.49 


Fe n 


2382.7652 


0.328 


10.41 



" Oscillator strengths taken from Morton (1991). 
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Table 5. Frequency distribution fitting parameters 



Form 


Grange 




logK 


a 




logNj 




Power-Law 


[2.55 


- 5.03] 


6.59±2.12 


-1.384d 


:0.105 








[2.55 


- 3.40] 


5.55±2.88 


-1.330d 


:0.144 








[3.40 


-3.83] 


4.03±2.01 


-1.056d 


:0.199 








[3.83 


- 5.03] 


4.38±2.39 


-1.279d 


:0.118 






Double Power-Law 


[2.55 


- 5.03] 


-21.78±0.31 


-1.162d 


:0.118 


20.60+0.24 


-2.07+0.51 




[2.55 


- 3.40] 


-22.64±1.25 


-1.327d 


:0.307 


21.19±1.44 


-5.98+2.07 




[3.40 


-3.83] 


-21.28±0.19 


-0.793d 


:0.233 


20.45±0.14 


-2.49+0.74 




[3.83 


- 5.03] 


-21.63±0.25 


-0.920d 


:0.199 


20.51±0.17 


-2.10+0.44 


Gamma 


[2.55 


- 5.03] 


-21.95±0.36 


-l.OlOd 


:0.172 


20.93+0.18 






[2.55 


- 3.40] 


-22.53±1.54 


-1.185d 


:0.329 


21.28+0.89 






[3.40 


-3.83] 


-20.96±0.28 


-0.485d 


:0.324 


20.45+0.16 






[3.83 


- 5.03] 


-21.80±0.35 


-0.853d 


:0.198 


20.82+0.18 





Table 6. Absorption Distance Path - Data for Figure [8] 



■Grange 




<z> 


A'qso 


A'dla 


A'subDLA 


AX 


nDLA(10-') 


f^DLA+subDLA(10') 


2.55- 


3.40 


3.168 


11 


12 


23 


125.922 


1.43 + 0.33 


1.71+0.33 


3.40- 


3.83 


3.618 


78 


17 


19 


125.922 


1.41+0.26 


1.65 + 0.26 


3.83- 


5.03 


4.048 


11 


11 


18 


125.922 


0.97 + 0.22 


1.21+0.22 
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Appendix : Notes on individual systems 

In this Appendix we discuss systems where we note differences between our measurements and measurements by others. 

1. PSS 0131+0633 (zem= 4.430). Peroux et al. (2001) report two sub-DLA candidates at Zabs= 3.17 and 3.61 with log NQi i) = 
19.9 and 19.8 respectively. For the second system however, Peroux et al. (2001) give in Table 5 a redshift of Zabs - 3.61 for 
the H I line and Zabs = 3.609 for the metal hues but a redshift of Zabs = 3.69 and log A'(H i) = 19.5 in their Section 8 namely, 
"Notes on Individual Objects". We also detect the first sub-DLA at Zabs = 3.173 with logA^(H i) = 19.95 + 0.10. Metal lines 
at that redshift are detected in the red part of the spectrum. For the second DLA candidate, we confirm the presence of the 
Zabs = 3.689 absorption with log A^(H i) = 19.50 ± 0.10 in agreement with the sub-DLA candidate reported in Section 8 of 
Peroux et al. (2001). Metal lines at that redshift are also observed in the red part of the spectrum. 

2. PSS0747+4434 (zem = 4.435). Peroux et al. (2001) report two DLA candidates at z^mahs = 3.76 and 4.02 with log A'(H i) = 
20.3 and 20.6 respectively, which are confirmed by our observations. We measure log A^(H i) = 20.00 ± 0.20 at Zabs = 3.762 
and log A'^(H i) = 21.10 ± 0.15 at Zabs = 4.019. No metal Unes at Zabs = 3.76 are observed in both surveys. We remind the 
reader that our data are of higher spectral resolution. 

3. SDSS 0756+4104 (zem = 5.09). We detect a sub-DLA candidate at Zabs = 4.360 with log N(U i) = 20.15 ± 0.10. Although 
no Lyman series and metals are detected for this system, the spectrum has a high enough SNR to fit the H i line well and to 
ascertain that this aborber is highly likely to be Damped. 

4. BR 0951-450 (zem = 4.35). PMSI03 report the detection of two DLA candidates at Zabs = 3.8580 and 4.2028 with log A'(H i) 
= 20.6 and 20.4 respectively, which are confirmed by our observations. We measure log A^(H i) = 20.70 ± 0.20 at Zabs = 3.856 
and log A'(H i) = 20.35 + 0.15 at Zabs = 4.202. We discover an additional Damped candidate at Zabs = 3.235 with log A^(H i) 
= 20.25 ± 0.10. We believe that this system is not included in the statistical sample of PMSI03, because it falls below the 
threshold of log A'^(H i) = 20.30. Metal Unes at that redshift are observed in the red part of the spectrum. 

5. BRl 1013+0035 (Zem = 4.38). PMS103 report one DLA candidate at Zabs = 3.103 with log A'(H i) = 21.1. This system is 
not included in our sample since its redshift is smaller than the minimum redshift beyond which a DLA/sub-DLA could be 
detected in our spectrum. One Damped candidate was discovered at Zabs = 3.738 with log A'^(H i) = 20.15 ± 0.20. We believe 
that this system is not included in the statistical sample of PMS103 because it falls below the threshold of log A^(H i) = 20.30. 
Although no lyman series and metals are detected for this system, the SNR in the spectrum is high enough to ascertain that 
this absorber is highly Ukely to be Damped. 

6. PSS 1057+4555 (Zem = 4.126). Peroux et al. (2001) report three DLA candidates at Zabs = 2.90, 3.05 and 3.32 with log A'(H i) 
= 20.1, 20.3 and 20.2 respectively, which are confirmed by our observations. We measure log A'(H i) = 20.05 ± 0.10 at 
Zabs = 2.909 , log NiU i) = 19.80 ± 0.15 at Zabs = 3.058 and log NQl i) = 20.15 ± 0.10 at Zabs = 3.317. Metal Unes at 
all three redshifts are observed in the red part of the spectrum. We report a new sub-DLA at Zabs = 3.164 with log A'^(H i) = 
19.50 ± 0.20. Several metal lines are detected for this sytem. 

7. PSS 1253-0228 (zem = 4.007). Peroux et al. (2001) report two DLA candidates at Zabs = 2.78 and 3.60 with log NiU i) 
= 21.4 and 19.7, respectively. The system at Zabs - 2.78 is not included in our sample because its redshift falls below the 
minimum redshift beyond which DLA/sub-DLA could be detected in our spectrum. We measure log A'(H i) = 19.95 + 0.10 
at Zabs = 3.608. Several metal lines are detected for this sytem. 

8. J 1325+1123 (Zem = 4.400). We have discovered two sub-DLA candidates at Zabs = 3.723 and 4.133 with log NiU i) = 
19.50 + 0.20 and 19.50 ± 0.20, respectively. C iv metal lines are detected at Zabs - 3.723. This two sub-DLAs can be found 
in the list of SDSS DR5 systems available on Jason Prochaska's website. 

9. PSS 1633+1411 (z^,„em = 4.349). Peroux et al. (2001) report one sub-DLA candidate at Zabs = 3.90 with log N(}i i) = 19.8. 
Here, we measure log A^(H i) = 19.0 + 0.20 at Zabs = 3.909. We report a new DLA at Zabs = 2.880 with log A^(H i) = 
20.30 ± 0.15. Several metal lines are detected for this sytem. 

10. PSS 1646+5514 (zem - 4.084). No DLA candidate is detected by Peroux et al. (2001) along this Une of sight. We report here 
two sub-DLA candidates at Zabs = 2.932 and 4.029 with log A^(H i) = 19.50 ± 0. 10 and 19.80 ± 0. 15, respectively. No metals 
are seen over the available wavelength coverage. 

11. PSS 2122-0014 (zem = 4.084). Peroux et al. (2001) report two candidates at Zabs = 3.20 and 4.00 with log N(n i) = 20.3 and 
20.1, respectively. Here, we measure log A'(H i) = 20.20 + 0.10 and 20.15 + 0.15 at Zabs = 3.207 and 4.001. Several metal 
Unes are detected for both systems. We report a new sub-DLA candidate at Zabs = 3.264 with log A'^(H i) = 19.90 ± 0.15. 
Several metal lines are detected also for this system. 

12. PSS 2154+0335 (zem = 4.359). Peroux et al. (2001) report two candidates at Zabs = 3.61 and 3.79 with log A^(H i) = 20.4 and 
19.70, respectively. No candidate with log A'^(H i) > 19.50 is detected in our spectrum. 

13. PSS 2155+1358 (Zem= 4.256). Peroux et al. (2001) report one DLA at Zabs = 3.32 with log NiU i) = 21.10. Dessauges et al. 
(2003) report three sub-DLAs at Zabs = 3.142, 3.565 and 4.212 with log A^(H i) = 19.94, 19.37 and 19.61. Here, we measure 
log NiU i) = 19.90 ± 0.20, 20.75 ± 0.20 and 20.00 ± 0.10 at Zrmabs = 3.143, 3.318 and 4.21 1, respectively. Several metal lines 
are detected for these systems. 

14. PSS2344+0342 (zem = 4.340). Peroux et al. (2001) report two DLAs at Zabs = 2.68 and 3.21 with log NiU i) = 21.10 and 
20.9. The first system is not included in our sample because its redshift falls below the minimum redshift beyond which a 
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DLA/sub-DLA could be detected in our spectrum. We measure log A^(H i) = 21.20 ± 0.10 at Zabs = 3.220. Dessauges et al. 
(2003) report one sub-DLA at Zabs = 3.882 with log NQi i) = 19.50. We measure log Ar(H i) = 19.80 ± 0.10 at Zabs = 3.884. 
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Fig. .1. Spectra of the 100 confirmed Lyman-a absorption with log A'CH i) > 19.5. In each panel the Voigt profile fit corresponding 
to the best A^(H i) value is overplotted. The Lyman series absorption lines and one characteristic associated metal absorption 
feature (when metals are present over the available wavelength coverage and not blended with another lines) are presented in 
additional sub-panels. 
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